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ABSTRACT 

We report the detection of extended warm ionized gas in two powerful high-redshift radio galaxies, NVSS J210626-314003 at z=2.10 and 
TXS 2353-003 at z=1.49, that does not appear to be associated with the radio jets. This is contrary to what would be expected from the alignment 
effect, a characteristic feature of distant, powerful radio galaxies at z>0.6. The gas also has smaller velocity gradients and line widths than most 
other high-z radio galaxies with similar data. Both galaxies are part of a systematic study of 50 high-redshift radio galaxies with SINFONI, and are 
the only two that are characterized by the presence of high surface-brightness gas not associated with the jet axis and by the absence of such gas 
aligned with the jet. Both galaxies are spatially resolved with ISAAC broadband imaging covering the rest-frame R band, and have extended wings 
that cannot be attributed to line contamination. We argue that the gas and stellar properties of these galaxies are more akin to gas-rich brightest 
cluster galaxies in cool-core clusters than the general population of high-redshift radio galaxies at z<;2. In support of this interpretation, one of 
our sources, TXS 2353-003, for which we have Ha narrowband imaging, is associated with an overdensity of candidate Ha emitters by a factor 
of ~8 relative to the field at z=1.5. We discuss possible scenarios of the evolutionary state of these galaxies and the nature of their emission line 
gas within the context of cyclical AGN feedback. 
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1. Introduction 


Powerful radio galaxies at redshifts z>l-2 (HzRGs^ are par- 
ticularly massive (seve ral 10 11 M 0 , e.g.. ISevmour et ai1l2007t 


of 1000 Mr. : . yr 1 Archibald et al 

2001 

Reuland et al 

2004 

Seymour et al. 2012F" Ivison et al. 

2012 

; Barthel et al. 

2012 

Drouartet al. 2014|), relatively evolved (Nesvadba et al. 

2011a 


galaxies. Many are surrounded by overdensities of satellite 
galaxies, making them interesting tracers of overdense regions 


rich galaxy clusters (e.g.. Chambers et al. 

1996; Le Fevre et al. 

1996 

; Kurk et al. 2004a; Venemans et al. 

2007c Galametz et al. 

2012 

; Wvlezalek et al. 2013). Their co-moving number den- 


sity corresponds t o that of massive clusters, when corrected fo r 
duty-cycle effects (IVenemans et al.ll2007l : INesvadba et al.ll2008l) . 
Particularly detailed studies of individual HzRGs suggest that 
at least one subset of this population could be the progenitors 
of the central galaxies in massive clusters in the nearby Universe 
(e.g. jHatch et all2009t) . They follow the low-redshift Kormendy 
relation of powerful radio and brightest cluster galaxies when 
allowing for passive luminosity evolution ( Targett et al . 20LI). 
Many HzRGs are themselves gas rich (e.gJvan Oiik et™ 1 1 9971 


INesvadba et aljf2008ll0gle et al.l l201 2: Em onts et all2014h . and 

while they have no well established, virialized halos of hot gas 
(lOverzier et al.ll2005l) . HzRGs do hav e extended (few 100 kpc) 
reservoirs of diffuse warm ionized gas ( Villar-Martm_e^alj2003b 


with embedde d clouds or filaments of neutral (van Oiik et al 


Il997t lAdams et alJ l2009h and molecular gas dNesvadba et al 


* Based on observations carried out with the Very Large Telescope 
of ESO under Program IDs 079.A-0617, 084.A-0324, 085.A-0897, and 
090.A-0614. 


1 20091 iDe Breuck et alJ 2003: Emonts et al.ll2014h . As suggested 
by, e.g., Villar-Martm et al. ( 20031) these could be the vestiges of 
the gaseous reservoirs from which the brightest cluster and other 
galaxies formed. 


The central, high surface-brightness regions of HzRGs are 
well fitted with de Vaucouleurs profiles with relatively large 
radii, R~8 kpc dTargett et al.l 1201 fl) . but extended, often irreg¬ 
ular line and continuum emission in HzRGs_is found prefer¬ 
entially along the radio jet axis (e.g., ICimatti et al.li 19931) . The 
reasons for this alignment effect are still not fully understood. 
Polarimetry suggests that light scattering on dust grains could 
be the origin of at least some of the extended continuum light; 
other authors favor jet-tri ggered star formation dChambers et al.l 
1988: iKfamer et al.l 12005 ) as the cause of the extended contin¬ 
uum emission. 


The kinematics in the exten ded gas around HzRGs 
are often strongly pertur bed (e.g., iBaum & McCarthy! [2000; 

1 Villar-Martm et al. 12003 !). SINFONI imaging spectroscopy 
showed this is the case in individual regions of these galaxies, 
and through out their extended emission-line gas ( Nesvadba et al.l 
2006, 120081) . Warm ionized gas seems to be a major compo¬ 
nent of the ISM in these galaxies , with ionized gas masses o f 
up to > 10 1() M 0 in some cases (iNesvadba et al.l 2008 1 2006 ). 
comparab le to the molecular gas masses found in some HzRGs 

(e.gjEmonts et alJl2014t INesvadba et alJl2009l: IDe Breuck et al.l 

12005 1. Broad line widths and abrupt velocity offsets of up to 
2000 km s 1 found in this gas are not reconcilable with grav¬ 
itational motion. Instead they suggest that a small part of the 
kinetic energy release from the AGN, through powerful radio 
jets, accelerates large fractions of the ISM in these galaxies to 
velocities well above the escape velocity from the underlying 
dark-matter halo, thereby driving winds and removing the fuel 
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for subsequent star formation. This scenario is broadly supported 
by hydrodynamic models showing that jets may indeed acceler¬ 
ate relatively dense gas to the observed velocities of a few 100 
km s _1 dWagner et al.ll2012h . 

Using VLT near-infrared imaging spectroscopy we have car¬ 
ried out a large survey of the rest-frame optical line emission 
of the warm ionized gas in a set of 50 radio-selected HzRGs at 
z<;2 spanning three orders of magnitude in radio power from 
a few times 10 26 to a few times 10 29 W Hz 1 at 1.4 GHz in 
the rest frame, and more than two orders of magnitude in radio 
size, from 1 kpc to about 300 kpc (Collet et al. 2015, submit¬ 
ted, Nesvadba et al. 2015, in prep.). Emission-line regions are 
well aligned with the radio jet axis in all but two galaxies with 
extended line emission, which in addition to kinematic, energy, 
and timing arguments suggests that this gas represents outflows 
of warm ionized gas driven by the overpressurized cocoon of 
hot plasma inflated by the radio source. These galaxies a re de- 
scribed in greater detail in iNesvadba et al.1 (l2006l 120071 [2008, 
1201 lah and Collet et al. (2015). 

Here we discuss the two radio galaxies, TXS 2353—003 and 
NVSS J210626—314003, that do not follow these global trends. 
Their extended regions of warm ionized gas with similar surface 
brightness to the other sources show strong offsets in position 
angle to the radio jet axis, contrary to what is expected from 
the alignment effect and observed in the other galaxies within 
the SINFONI sample. Gas associated with the radio jet axis, if 
present at all, is not a distinctive component compared to gas 
along other directions in our SINFONI data. Gas that is seen 
in projection to be associated with the jet axis may or may not 
imply a physical association; however, it is clear that gas that 
is globally and significantly offset from the jet axis as projected 
on the sky is not reconcilable with a scenario where the gas is 
embedded in an overpressured cocoon inflated by the on-going 
jet activity. In addition, the line widths are low compared to the 
rest of our sample. Line widths near the nucleus are higher, and 
may be affected by the central radio source. 

TXS 2353-003 and NVSS J210626-314003 are among 
the sources with the largest radio sizes in our sample, with 
largest angular sizes LAS=38.8" and 24.2" at 4.86 GHz 
dDe Breuck et al.l 1200 1) and 2.368 GHz, respectively (corre¬ 
sponding to 310 kpc and 190 kpc at their redshifts of z=1.49 
and z=2.10, respectively). Both have very regular centimeter 
radio morphologies (Collet et al. 2015, Nesvadba et al. 2015 
in prep.), and do not have bright radio cores. Both are associ¬ 
ated with overdensities seen in projection along nearby lines of 
jiight of distant (z>1.3) IRAC sources in the CARLA survey of 
IWvlezalek et al.l (2013 ). TXS 2353-003 is within the densest re¬ 
gion of their overall sample of 200 HzRGs. We highlight the 
similarities between our galaxies and brightest cluster galaxies 
at low redshift, including the extended emission-line gas, which 
resembles the warm ionized gas in and around the central galax¬ 
ies of cool-core clusters in several respects. 

We describe our SINFONI and ISAAC observations and data 
reduction in i)2] and discuss the results of these observations in 
fj3]for the line emission and in fj4]for the continuum, where we 
also fit surface brightness profiles. In : [5] we present our narrow- 
band search for line emitters around TXS 2353-003 and argue 
that this galaxy is surrounded by an overdensity of line emitters, 
like many other HzRGs at somewhat greater redshifts. We com¬ 
pare these results with the IRAC results obtained as part of the 
CARLA survey C Wvlezalek et a l. 2013 ) in : [6] In : {7] we discuss 
several hypotheses regarding the nature of the emission-line re¬ 
gions in these two galaxies, which, given their strong misalign¬ 
ment relative to the jet axis, are unlikely to be directly associ- 



Fig. 1 . Results of our observations of TXS 2353-003. Top left: 
5.5 GHz contours from ATCA. Clockwise from top right to 
bottom left: SINFONI maps of [OIII]A5007 surface brightness, 
FWHM line widths, and relative velocities. The square in the top 
left panel shows the size of the SINFONI maps compared to ra¬ 
dio size. Black contours show the 5.5 GHz radio continuum, and 
the dotted square in the top right panel illustrates the region from 
which the bottom spectrum in Fig. [3] was extracted. The arrow 
in the lower left panel shows the direction of the radio jet, and 
the number gives the jet size in arcsec. Ellipses in the lower left 
corner of each panel show the FWHM size of the point spread 
function. 


ated with an expanding cocoon driving an outflow. We summa¬ 
rize our results in Sj8] Throughout our analysis we adopt a flat 
cosmology with Ho = 70 km s^ 1 Mpc -1 , £2a=0.7, Q.m = 0.3. 
For NVSS J210626-314003 at z=2.1, this corresponds to D/. = 
16.6 Gpc, and for TXS 2353-003 at z=1.49 to D L =10.8 Gpc, 
respectively; 1" corresponds to 8.4 kpc for both galaxies. 

2. Observations and data reduction 

2.1. VLT/SINFONI imaging spectroscopy 

We observed both galaxies with the near-infrared imaging spec¬ 
trograph SINFONI on UT 4 of the Very Large Telescope (VLT) 
of ESO. The NVSS J210626-314003 data were taken on 2009 
October 12, and those of TXS 2353-003 on 2010 August 13. The 
data were obtained as part of a comprehensive observational pro¬ 
gram to study the emission-line gas kinematics of 50 powerful 
radio galaxies at z>2. We used the H+K band covering a spec¬ 
tral range of 1.45-2.45 jim at a spectral resolving power R=1500. 
The data were taken during good to moderate observing condi¬ 
tions, with a typical point spread function of 1" x 1" along right 
ascension and declination, dominated by the seeing. We used the 
250 mas pixel scales with a field of view of 8"x8", covering a 
field of 67 kpc x 67 kpc at z=2. 

Total integration times were 145 min for TXS 2353-003 and 
120 min for NVSS J210626-314003, and were obtained in one- 
hour sequences of individual 300 s exposures. Our data reduc¬ 
tion relied on a combination of IRAF and IDL routines (for de¬ 
tails see, e.g., Nesvadba et al.ll20T la b). The absolute flux scale 
was measured by comparing our data with observations of stan- 
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Fig. 2. Results of our observations of NVSS J210626—314003. 
Top left: 5.5 GHz contours from ATCA. Clockwise from top right 
to bottom left: SINFONI maps of [OIII]A5007 surface bright¬ 
ness, FWHM line widths, and relative velocities. The square in 
the top left shows the size of the SINFONI maps compared to ra¬ 
dio size. The black contour marks the continuum position in our 
SINFONI cube. Ellipses in the lower left corner of each panel 
show the FWHM size of the point spread function. 


dard stars with known magnitudes from the 2MASS survey that 
were observed at the end of each observing session. We also de¬ 
termined the size of the seeing disk from these stars. 

2.2. VLT/ISAAC imaging 

Both galaxies were also observed with ISAAC on UT 3 of the 
VLT between 2012 September 30 and 2012 October 2 through 
the service-mode program 090.A-0614. TXS 2353—003 was ob¬ 
served in the H band and through the narrowband Fell filter 
(NB 1.64) with FWHM= 0.025/im (corresponding to a veloc¬ 
ity range of Av =4570 km s -1 ), so we cover Ha and [Nil] at 
z=1.49. NVSS J210626-314003 was only observed in the Ks 
band. Both broadband images cover roughly the rest-frame R 
band. 

With the NB 1.64 filter we obtained a total of 248 minutes of 
on-source observing time of TXS 2353—003 split into 99 ex¬ 
posures with individual observing times of 150 seconds, and 
grouped into observing sequences of one hour each. After visual 
inspection we discarded 17 of these frames because of residuals 
from a strongly saturated star. The on-source observing time was 
thus 205 minutes. 

In the H band we obtained 29 exposures for TXS 2353- 
003 with individual observing times of 72 seconds, correspond¬ 
ing to six detector readouts after 12 seconds. We discarded 
one frame which contained the trail of a satellite. In the K 
band for NVSS J210626—314003, we obtained 25 exposures 
of 90 seconds, corresponding to six detector readouts after 
15 seconds. Total observing times in the H and K bands were 
therefore 35 minutes and 38 minutes for TXS 2353—003 and 
NVSS J210626-314003, respectively. 

Data were dark-subtracted and flat-fielded, where we con¬ 
structed the sky flats directly from the science data. All frames 
were then combined with the IRAF task xdimsum. Flux cali¬ 


bration relied on standard stars observed during the same night 
and at similar air mass to the science data. We used the Starlink 
GAIA software ( [Draper et al.ll2014h to calculate the astrometry 
relative to 2MASS reaching rms=0.07" for TXS 2353-003, and 
rms=0.3" for NVSS J210626-314003, respectively. We also 
used 2MASS to verify our flux calibration, finding a good agree¬ 
ment within < 10% . The sizes of stars in the image suggest the 
seeing was FWHM=0.5 ;/ in the H-band and NB 1.64 image of 
TXS 2353-003, and of FWHM=0.4" in the Ks-band image of 
NVSS J210626-314003. The limiting rms surface brightness in 
these images is 26.2 mag arcsec -2 and 23.6 mag arcsec -2 in 
the narrow and broadband image of TXS 2353—003, respec¬ 
tively. For the Ks-band image of NVSS J210626-314003, we 
find 23.7 mag arcsec -2 . 

2.3. Narrowband observations of TXS 2353—003 

To construct a continuum-subtracted line image from the nar¬ 
rowband image of TXS 2353-003, we approximate the underly¬ 
ing continuum in the narrowband filter by a scaled version of the 
flux measured through the broadband filter, taking into account 
the relative filter bandwidths of the narrow and broadband fil¬ 
ters, and their transmission. Here we assume that the continuum 
is flat and uniform across the H-band filter. 

We measured the zero points in the narrowband filter from 
standard stars. From the zero points in each filter, filter band- 
widths, and transmissions of the narrow and broadband fil¬ 
ters, we determine that the continuum flux density observed 
through the narrowband image corresponds to about 7.2% of the 
flux density measured through the broadband image. We would 
therefore need to scale the broadband image by this value before 
subtracting it from the narrowband image to obtain the flux den¬ 
sity from the emission lines. When directly comparing the flux 
density of 12 stars in both images, we empirically find that a very 
similar scaling factor of 7.6% is best to minimize the residuals of 
the stellar continuum after the subtraction, and this is the factor 
that we used. 

The redshifted wavelength of Ha in the radio galaxy is 
slightly blueshifted relative to the wavelength of peak through¬ 
put of the narrowband filter. Comparing this with the flux mea¬ 
sured in the SINFONI data cube, and assuming that the rel¬ 
ative flux calibration between both images is perfect, we find 
that we are missing between 20 and 25% of the line flux in the 
continuum-subtracted line image. For other Ha emitters ([J5]), for 
which we have no precise redshift estimates, we would first need 
to obtain spectroscopy before we can make similar corrections. 

To evaluate the completeness limit of our data, we added 
faint artificial sources with the size of the seeing disk at sky 
positions in each frame of the narrow- and broadband images 
before re-reducing these data and extracting sources in the same 
way as for the final scientific analysis. We used fluxes between 2 
and 11 times the rms in the central regions of the image. Adding 
the artificial sources to the individual frames rather than the final 
image helps to take misalignments between frames into account. 
We find that we recover 90% of sources with surface bright¬ 
ness above 1.5 x 10 -19 erg s -1 cm -2 A -1 arcsec -2 (or above 
23.6 magAB arcsec -2 ). 


3. Imaging spectroscopy of the warm ionized gas 

The SINFONI maps of NVSS J210626-314003 and 
TXS 2353-003 are shown in Figs. Q] and [2] respectively. In 
NVSS J210626—314003 at z=2.104, [OIII]AA4959,5007 and 
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HJ3 fall into the H band, and Ha, [NII]AA6548,6583, and 
[SII]AA6716,6731 fall into the K band, both of which we cover 
with our data cubes. In TXS 2353—003, we identify Ha and 
[NII]AA6548,6583 at wavelengths corresponding to a redshift 
of z=1.49. At this redshift, [OIII]AA4959,5007 and H/3 fall 
at redshifts below the lower wavelength cutoff of our grating, 
Amin = 1 -45 gm. 

Our redshift measurement for TXS 2353-003 does not agree 
with the previous redshift of z=2.59 measured with longslit 
spectroscopy, which was thought t o have detected Lyq i n the 
observed optical wavelength range (iDe Breuck et al.ll200ll) . but 
was in hindsight probably affected by a Cosmic ray. We did not 
identify any possible emission line consistent with the previous 
redshift z=2.59, although [OIII]AA4959,5007 and H/3 should 
fall into the H band, and [NII]AA6548,6583, and Ha should fall 
into the K band at that redshift. The rest-frame UV spectrum 
shows no other line candidate consistent with z=2.59; the value 
for Lya at the redshift z. = 1.49 we measure here falls well below 
the UV atmospheric cutoff. 

Both galaxies have bright, spatially well-resolved line emis¬ 
sion. We fitted the extended line emission in both galaxies with 
single Gaussian profiles. Spectra were extracted from small aper¬ 
tures of 0.4" x 0.4", which helps to reduce the strongest pixel-to- 
pixel noise while safely oversampling the seeing disk so that no 
spatial information is lost. The data were also smoothed in the 
spectral direction by three pixels, without loss of spectral res¬ 
olution. Our fitting routine interpolates over wavelength ranges 
dominated by bright night-sky line residuals, and we only map 
spatial pixels where a line was detected at a S/N >3. 

Line emission extends beyond the bright continuum seen 
in the SINFONI cubes (but not in the ISAAC broadband im¬ 
age which we will discuss later in more detail and show 
in Figs. □ and [8] which is deeper in the continuum). The 
sizes of the extended emission-line regions of NVSS J210626- 
314003 and TXS 2353-003 are 5.3"x2.0" (44 kpcx 17 kpc) and 
2.5"x2.0"(21 kpcx 17 kpc) along the major and minor axis, re¬ 
spectively, and down to surface-brightness limits of 4.3 x 10 17 
erg s"" 1 cm" 2 arcseW 2 and 6.2 x 10“ 17 erg s"" 1 cm" 2 arcseW 2 . 

In NVSS J210626—314003 the emission-line region is 
roughly centered on the continuum peak. In TXS 2353—003, the 
gas appears more lopsided. This is also seen in the continuum- 
subtracted ISAAC narrowband image (source a in Fig. lA. Ik The 
morphology of the ionized gas in both observations is fairly con¬ 
sistent: a bright emission-line region associated with the radio 
galaxy, more extended toward the southeast and northwest. The 
integrated spectrum extracted from this fainter region, shown as 
a dashed region in Fig. Q] clearly shows Ha and [Nil] A6583 
(Fig. 01 . Fluxes and other line properties are listed for both galax¬ 
ies in Table [I] In TXS 2353—003, we note a faint continuum 
emitter along the direction of extended gas, but at a larger dis¬ 
tance. In NVSS J210626—314003 no further continuum source 
is observed along the major axis of the emission-line region. 

In the SINFONI cube, NVSS 210626-3134003 shows a 
smooth velocity gradient of about 600 km s 1 over 25 kpc. In 
TXS 2353—003 we observe a more irregular velocity pattern, 
which could be at least in part due to the lower signal-to-noise ra¬ 
tio and irregular line profiles. In NVSS J210626-3134003, where 
we observe [OIII]A5007 and Ha at good signal-to-noise ratios, 
the gas kinematics and morphologies are similar in both lines. 

Both galaxies have broad line widths with 
FWHM~800 km s" 1 around the continuum peak which 
we associate with the central regions of the galaxies. These lines 
are broader than those of star-forming galaxies at similar red- 
shifts, including the most massive, intensely star-forming dusty 



Wavelength (A) 



Wavelength (A) 


Fig. 3. Top: Spectrum of TXS 2353—003, integrated over all spa¬ 
tial pixels where the Ha emission from the galaxy is detected 
at > 3(7, and corrected for velocity shifts. Bottom: Integrated 
spectrum of the northeastern part of the emission-line region of 
TXS 2353—003, defined by the dotted square in Fig.Q] 


galaxies such as submillimeter galaxies (e.g., ISwinbank et al.l 
120061) . and are consistent with those found in other powerful 
radio galaxies at similar redshifts (e.g.. iNesvadba etakl f2008: 
ICollet et akl 120151) . Following the arguments given in these 
studies, we consider the broad line widths as evidence that the 
gas in the central regions of our sources is stirred up by the 
transfer of kinetic energy from the radio source. 

The gas kinematics at larger radii are more quiescent in both 
sources. Line widths in the extended gas of NVSS J210626- 
314003 and TXS 2353-003 are FWHM-300-400 km s“' 
(Figs. Q] and [2}, more akin to those observed in a quasar illu¬ 
mination cone at z~2 dLehnert et a l. 2009) and Lyq blobs with¬ 
out a powerful cent ral radio source (e.g.. IWilman et aU 12005: 
lOverzier et al.ll2013l) . than to powerful radio galaxies, including 
the central regions of our two sources here. 

A marked difference to other HzRGs with similar data is that 
most of the emitting gas is at large position angles from the ra¬ 
dio jet axis. Typically, the jet axis and the semi-major axis of 
the gas are aligned within about 20°-30° (Fig. [4] see also Collet 
et al. 2015, Nesvadba et al. 2015 in prep.). Given the size of 
the minor axis of the extended line emission of up to ~ 15—20 
kpc and beam smearing effects in the radio and near-infrared, 
this range is on the order of what would be expected if gas and 
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PA jet 

Fig. 4. Position angle (north through east) of the gas (ordinate) 
as a function of the position angle of the radio jet for our two 
sources. The outer boundaries of the red, blue, and green stripes 
indicate offsets of 20°, 45°, and 90°, respectively. Red stars in¬ 
dicate our two sources discussed here, small black dots show 
the remaining HzRGs with SINFONI data and well-determined 
position angles (Nesvadba et al. 2015, in prep.). 


jet were associated with each other, for example, if the jet in¬ 
flated a hot bubble of overpress urized gas like in the cocoon 
model JBegelman & Cioffil 1989 ). In NVSS J210626-314003 
and TXS 2353—003; however, the relative position angle be¬ 
tween extended gas and radio jet axis is much larger, 60° and 
90°, respectively, and the jet is not embedded within the emis¬ 
sion line gas, as expected from a cocoon inflated by the jet. 
Figure[4]shows that both galaxies clearly stand out from the over¬ 
all sample of HzRGs with SINFONI observations. 


3.1. Line diagnostics 


The ratios of bright optical line emission provide important con¬ 
straints on the ionizing source and other gas p roperties in the 
extended emission-line regions of galaxies (Baldw in et al. 1981; 
Veill eux & Osterbrockl 119871 : IKewlev et all 2006). Collet et aU 
( 20151) discussed the emission-line ratios of NVSS J2106- 
314003 in the context of their overall sample of moderately pow¬ 
erful HzRGs, finding [NII]/Ha and [OIII]/H/3 ratios that cor¬ 
respond to neither the classical starburst nor the AGN branch. 
They investigate several potential reasons for this, including low 
metallicites, relatively diffuse gas clouds, high gas pressures, 
or a mixture of heating from shocks and photoionization. In 
TXS 2353-003, we only have [Nil] and Ha measured with a 
fairly high ratio of [NII]/Ha=0.7 in the extended gas (Fig. [3]), 
which alone is sufficient to safely attribute the line emission to 
either shocks or photoionization from the AGN (Fig. 0. Off- 
axis gas not aligned with the radio jet axis and out to a galac- 
tocentric radius of 65 kpc in the radio galaxy 3C265 at z=0.8 
may indeed be powered by shocks, perhaps with additional pho¬ 
toionization from_soft X-ray emission from the precursor of the 


shock (So lorzano-Ifiarre a et al. 20021), However, as argued by 


iLe Tiran et alJ (1201 ll) and Nesvadba et al.l (i2006h . it is difficult 


to reach high Ha surface brightnesses as observed with shocks. 
The same is the case for shock-heated, likely infalling clouds in 
the cluster central galaxy NGC 4696 dFarage et alJl20ld) . 



log([NII]/H-alpha) 


Fig. 5. Diagnostic diagram showing the ratio of [OIII]/H/3 
as a function of [NII]/Hq line ratio dBaldwin et al.l 119811; 
[Veilleux & Osterbrockl 119871) . The red dot shows the position 
of NVSS J210626-314003. The red hatched area shows the 
region spanned by the [NII]/Ha ratio in the off-nucleus gas 
in TXS 2353-003. Light blue dots indicate the distribution of 
line ratios for galaxies in the local Universe (from SDSS; see 
iKauffmann etalj 2003), and t he dashed blue line shows the di¬ 
agnostic of lKewlev et al.l (120061) to separate HII regions and star¬ 
forming galaxies (below) from AGN (above). The solid blue 
lines indicate how the distribution of line ratios for both SF 
galaxies and AGN are expected to evolve out to high redshift 
(iKewlev et al.ll2013l their model 4 with metal-poor narrow-line 
regions): the two HzRGs are clearly consistent with the predic¬ 
tions for high redshift AGN. 


Like TXS 2353-003, NVSS J210626-314003 also falls out¬ 
side the star formation branch; however, it also falls outside the 
sequence formed by low-redshift AGN in the SDSS. The line 
ratios are consistent with low-metallicit y gas photoionized by 
an AGN, as expecte d in the models of iGroves et al] (120061) and 
IKewlev et al] (120131) . and as shown with the solid blue lines in 
Fig. E These lines show the range of line ratios expected for 
galaxies at z=l .5 with intense star formation typical for galaxies 
at these redshifts, and AGN with metal-poor narrow line regions 
and a metal enrichment history predicted by cosmological mod¬ 
els of galaxy evolution (Dave et al. 20H). These conditions are 
referred to as “model 4” by iKewlev et alJ (i2013l) . J-band spec¬ 
troscopy would be required to investigate whether TXS 2355- 
003 falls onto the low-redshift AGN branch or above. 


4. Continuum morphologies 

4.1. Line-free continuum images from SINFONI 

We also constructed line-free continuum images from our 
SINFONI cubes (shown as contours in Figs. [T] and [2]) for both 
galaxies by collapsing the data cubes along the spectral direc¬ 
tion over all wavelengths that are free from emission lines and 
prominent night-sky line residuals. We consider the continuum 
peak in the cube to be the location of the galaxy nucleus, which 
may be dominated either by stellar light or the AGN. At a pro¬ 
jected distance of about 3" (24 kpc) northeast of TXS 2353-003, 
we find another, fainter H-band source. We do not detect any 
line emission from this galaxy, so we cannot determine whether 
it is physically associated with the radio galaxy, or an inter¬ 
loper along the line of sight. NVSS J210626—314003 appears 
to be an isolated source, except for two much fainter sources at 
1.7" and 3.9" toward the southeast (corresponding to projected 
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separations of 14.3 kpc and 32.8 kpc at z=2.1, respectively), 
which however are not aligned with the extended emission-line 
regions. We do not detect any line emission from the nearer 
source, and the more distant source falls outside the field of view 
of SINFONI. It is therefore not clear if they are associated with 
the radio galaxy. For TXS 2353-003, for which we do not have 
broadband K-band imaging, we measure a K-band magnitude of 
18.1 ±0.2 mag from our SINFONI data cube within a 3" aper¬ 
ture. 


4.2. ISAAC broadband imaging 

In Figs. [7] and [8] we have already shown the rest-frame R-band 
morphologies of TXS 2353-003 and NVSS J210626-314003 
observed in the H and Ks band, respectively. Both galaxies 
are clearly spatially resolved, even with ground-based, seeing- 
limited images. This is best seen in Fig. [6] which shows the az- 
imuthally averaged surface-brightness profiles of both sources 
compared to the profiles of nearby stars taken from the same im¬ 
ages, showing a clear excess of continuum emission out to large 
radii in both galaxies. 

Both galaxies have a bright central region (which is more 
regular in NVSS J210626-314003 than in TXS 2353-003) 
surrounded by extended, low-surface-brightness halos. This 
is remarkable, given that the general population of passively 
evolving early-type galaxies at high redshift tends_to be more 
compact than their local analogs (e.g. , Daddi et all f2 005: 


Dokkum et al. 1120081: ISchawinski et al.ll201 1; Whitaker et al 


20121) . Furthermore, although about 1/3 of HzRGs are still ac¬ 
tively star-forming dDrouart et al.1 12014 ), their molecul ar gas 


reservoirs of a few 10 lu M 0 (e.g., Emonts et al. 2014) com¬ 
pared to few 10 11 M m of stellar mass ( Seymour et al J f2007 : 
iDe Breuck et al.1 1201 Oh suggest that most of their stella r mass 
growth has already been completed. [jargett et al] (1201 ll) found 
extended broadband emission when stacking 13 HzRGs at red- 
shifts between z=1.5 and 2.0 (as did Best et alii 19981 at slightly 
lower redshifts). However, only lHatch et al.1 (12009 ) have so far 
reported a qualitatively similar surface-brightness profile in an 
individual radio galaxy, MRC 1138—262 at z=2.16. Previous 
studies, e.g. JPentericci et al, (2001), have drawn attention to the 
irregularity of broadband morphologies in many cases, which 
they attributed to the presence of on-going mergers. 


4.3. Surface brightness profiles 

Both galaxies are spatially resolved into four or five resolution 
elements in the ISAAC images, which enables us to do a basic 
analysis of their surface-brightness profiles shown in Fig. [6] The 
goal of this analysis is to demonstrate that the extended low- 
surface brightness regions seen in the broadband images are not 
just the faint extensions of the surface-brightness profiles of the 
central regions, but separate components of the light profiles. 
We are mainly interested in whether (and over what radii) the 
profiles are consistent with the r 1 / 4 or de Vaucouleurs profile 
typical of early- type galaxies , corresponding to a Sersic profile 
with index n=4 (ISersicll 1963l 119681) . 

The profiles shown in Fig. [6] were obtained with the IRAF 
task Ellipse (IFreudlingll 19931) . measuring the surface bright¬ 
ness in circumnuclear elliptical annuli with radii increasing in 
steps of 0.1". Ellipse fits the ellipticity of these annuli, find¬ 
ing between 0.05 and 0.4 in our case. Solid red lines in Fig. [6] 
show the expected curves of n=4 Sersic profiles that match the 
central 1" of the light profile, with effective radii R=261 kpc 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 

R 1,4 [kpc] 



Fig. 6. Surface-brightness profiles of NVSS J210626-314003 
(top) and TXS 2353—003 (bottom) extracted from ISAAC H- 
band and Ks-band images, respectively. Red dots show the mea¬ 
sured azimuthally averaged surface brightnesses as a function 
of R 1 4 . Red lines show a Sersic fit with n —4, corresponding to 
the de Vaucouleurs law for classical elliptical galaxies, and con¬ 
volved with the point spread function in each image. The dashed 
black lines show the size of the seeing disk. Dotted horizonal 
lines show the 3 o limiting surface brightness of our images. 
Excess flux at large radii compared to the R 1 / 4 law is charac¬ 
teristic for brightest cluster galaxies in evolved galaxy clusters 
(e.g.. Schombe rtlll98'7l) . 


and R=15 kpc for TXS 2353-003 and NVSS J21062-314003, 
respectively. These effective radii are consistent with those we 
found with two-dimensional fits using Galf it dPeng et al.l2002i 
l2010i) . which we ran on images that were slightly smoothed to 
minimize pixel-to-pixel noise (by convolution with a 3 x 3 pixel 
Gaussian) and where we left all parameters free, except impos¬ 
ing that n=4. 

The comparison between the expected profiles for n=4 and 
the observations shown in Fig. [6] illustrates that only the in¬ 
ner regions of NVSS J210626—314003 can be fitted with a 
de Vaucouleurs profile, with radius R e = 15 kpc. However, 
the extended halo produces residuals greater than 3cr. For 
TXS 2353—003, the surface-brightness contrast between the 
central regions and the extended outer halo is lower, result¬ 
ing in a very large size of R e = 261 kpc for a de Vaucouleurs 
profile, which matches the central and outermost isophotes, 
but also leads to significant residuals at intermediate radii. 
This is in agreement with the visual impression that only 
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Fig. 7. left: ISAAC H-band broadband image of TXS 2353-003 
convolved with a two-dimensional Gaussian to match the spa¬ 
tial resolution of the Ha + [Nil] map shown in the right panel. 
right: Ha morphology of TXS 2353—003, convolved with a 
two-dimensional Gaussian to emphasize the morphology of the 
faint northeastern extension at a spatial resolution of 1.2". The 
color bar shows the contamination to the flux density measured 
in the broadband filter from Ha line emission in units of 10 19 
erg s’ 1 cm -2 A -1 arcsec 1 . Contours show the broadband K- 
band morphology and are identical to those shown in the left 
panel. Starting from the outermost contour, the broadband fluxes 
represented by each contour are 1.1,1.8, 2.5, and 5.4 x 10 -19 erg 
s -1 cm -2 A -1 arcsec -1 (3, 5, 7, and 15a in the convolved im¬ 
age, where a = 3.9 x 10 -2 ° erg s -1 cm -2 A -1 arcsec -1 ). The 
circles in the lower left of each panel show the size of the con¬ 
volved PSF. 


NVSS J210626—314003 has a distinctive, high-surface bright¬ 
ness core, whereas TXS 2353—003 is more diffuse. High- 
resolution imaging with the Hubble Space Telescope, for exam¬ 
ple, would be needed for a full analysis of the surface-brightness 
profiles in the inner regions of our galaxies. Nonetheless, these 
results already show that both galaxies show important depar¬ 
tures from simple de Vaucouleurs profiles at intermediate and 
large radii. Such extended wings are characteristic of central 
cluster galaxies (e.g.. lSchombertIl987lf . 

4.4. Line contamination 

Line emission in high-z radio galaxies can reach high equiv¬ 
alent widths and might therefore be a significant contam¬ 
inant in morphology measurements through b r oadband fil¬ 
ters (jPentericci et afll2001t iNesvadba et al.ll2008t iTargett et~ail 
1201 ill . We will now investigate the significance of this contami¬ 
nation in our analysis. 

(■.From our SINFONI imaging spectroscopy we know the 
line fluxes of the most prominent optical emission lines in both 
galaxies, Ha, [NII]AA6548,6583, and [OIII]AA4959,5007. In 
TXS 2353-003, the combined flux of all emission lines in the 
H band reaches an average surface brightness of 1.6 x 10 -15 
erg s -1 cm -2 arcsec -2 in a 2.15" aperture around the center 
of the galaxy, and 3 x 10 -16 erg s -1 cm -2 arcsec -2 along the 
northeastern extension. Outside this region, line emission is not 
detected and must therefore be below the 3 a upper surface- 
brightness limit of 1 x 10 -16 erg s -1 cm -2 arcsec -2 . 

The broadband morphology of TXS 2353—003 is shown in 
the left panel of Fig. |7] as a gray scale image and contours. The 
image was convolved with a two-dimensional Gaussian to a res¬ 
olution of 1.2", to match the size of the seeing disk of the Ha 
and [NII]A6583 line image shown in the right panel of Fig. [7] 




Fig. 8. left ISAAC K-band broadband image of NVSS J210626- 
314003 smoothed to the same spatial resolution as the 
SINFONI imaging spectroscopy, 1". right Ha morphology of 
NVSS J210626-314003. The color bar shows the contamination 
to the flux density measured in the broadband filter from Ha line 
emission in units of 10 -19 erg s -1 cm -2 A -1 arcsec -1 . Contours 
show the broadband K-band morphology and are the same as in 
the left panel. The contour levels correspond to 3, 5, 7, 15, and 
30 x the root-mean square of the broadband image smoothed to 
the SINFONI resolution of 1", which has r.m.s.= 3 x 10 -20 erg 
s -1 cm -2 A -1 arcsec -1 . The circles in the lower left of each 
panel show the size of the seeing disk. 


and to enhance the contrast of the faint extended emission in the 
broadband image against the background noise. 

The right panel of Fig. |7] shows the Ha+[NII] line image 
of TXS 2353—003, obtained by collapsing the data cube along 
the spectral direction over wavelengths where Ha+[NII] is de¬ 
tected. This makes the line morphology measurement more ro¬ 
bust than the maps obtained from Gaussian fits to each individual 
pixel, in particular in the faint outer parts of the emission-line re¬ 
gions, and allows the measurement of the upper limits from the 
same image that is also used to estimate the emission-line sur¬ 
face brightness. Dividing by the 2700 A width of the H-band fil¬ 
ter of IS AACQ we find that line emission contributes on average 
2.75 x 10 19 erg s 'em 2 arcsec 2 A 1 to the measured flux 
density in the H band near the center of TXS 2353—003, and 
of 1.1 x 10 -19 erg s -1 cm -2 arcsec -2 A -1 in the northeastern 
periphery. Outside this extended emission-line region, the flux 
densities from emission lines drop to below 3.7 x 10 20 erg s 1 
cm -2 arcsec -2 A -1 (our 3a limit). H-band surface brightnesses 
at the location of the extended Ha emission-line region are 
above 2.5 x 10 -19 erg s -1 cm -2 A -1 arcsec -2 , which suggests 
that line contamination does not dominate the overall broadband 
morphology of TXS 2353-003. 

The broadband image of NVSS J210626—314004 is shown 
in Fig. [S] as a gray scale image with contours (left panel), and 
as contours on top of the Ha and [Nil] line image (right panel), 
which was constructed in the same way as described above for 
TXS 2353—003. The line emission is closely approximated by 
Gaussian line profiles, and we therefore use the emission-line 
map derived from those fits and shown in Fig. [2] In a 2.15" 
aperture around the nucleus, we find an average surface bright¬ 
ness of Ha and [Nil] of 6.9 x 10 16 erg s 1 cm 2 arsec 2 and 
2.2 x 10 -16 erg s -1 cm -2 arcsec -2 in the periphery, respectively. 
This corresponds to a surface flux density of 2.3 x 10 19 and 
9.4 x 10 -20 erg s -1 cm -2 A -1 arcsec -2 , respectively, after di¬ 
viding by the filter width of the ISAAC K-band filter, 3000 A 1 . 


1 http://www.eso.org/sci/facilities/paranal/decommissioned/ 

isaac/inst/isaac_img.html 
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Outside the ridge of extended line emission, the K-band flux den¬ 
sities from the lines drop to below 3 x 10 20 erg s" 1 cm -2 A -1 
arc sec 2 (our 3cr limit on the emission-line surface brightness). 

Line emission in NVSS J210626—314003 extends over re¬ 
gions of broadband isophotes of 1.5 x 10“ 19 erg s^ 1 cm~ 2 and 
brighter (Fig. This comparison suggests that, outside the 
brightest areas seen in the Ha line image near the galaxy cen¬ 
ter, line emission is not a major contaminant of the broadband 
morphology in NVSS J210626-314003. 

Th e di ffuse Lyq halos discove red by IVillar-Martm et all 
(l2002ll and lVillar-Martm et al.1 (120031) would also produce fainter 
emission-line surface brightnesses than we observe here. For 
a typical SBi ya = 2 — 3 x 10~ 17 erg s -1 cm“ 2 arcsec~ 2 , we 
would expect Ha surface brightnesses between 2 — 5 x 10 1 x 
erg s - 1 cm 2 arcsec 2 , assum ing typical Lya/Ha ratios of 
8—13 ( Villar-Martin et al.ll2003b . Dividing by the 2700A filter 
width of the ISAAC H-band filter, this would correspond to 
0.7 — 1.5 x 10 -21 erg s -1 cm~ 2 A" 1 arcsec^ 2 , much lower than 
the surface brightness seen in the broadband filter. 



Fig. 9. Broad- to narrowband color-magnitude diagram of our 
candidate Ha emitters. Ha candidates are labeled with col¬ 
ored symbols, depending on the significance of their detection. 
Letters refer to the sources shown in Fig. IA.ll Source “a” is the 
HzRG. 


5. Ha candidates surrounding TXS 2353-003 

Narrowband imaging is a standard way of identifying emission¬ 
line galaxies at high redshifts, and has been successfully 
used to identify over-densities of actively star-forming galax- 
ies around high-redsh i ft galaxies for almos t two decades (e.g. . 


Chambers et al. 

1996 

Le Fevre et a 

1. 199 

it Kurk et al. 2004a; 

Venemans et al. 

2007 

Hatch et al. 

2011; 

Kovama et al. 2013; 

Havashi et al. 2012; Cooke et al. 20141). Most of these observa- 


tions were carried out for HzRGs at z>2, where Ha falls in the 
K band, which is relatively free from telluric night-sky lines. 


We also obtained Ha narrow-band imaging of TXS 2353- 
003, where Ha falls fortuitously into the NB 1.64 filter (al¬ 
though not perfectly into the center). Our main goal was to 
search for possible bright emission-line gas outside the small 
SINFONI field of view of 8”x8”, that would be associated with 
the environment rather than the radio galaxy itself. This also en¬ 
ables us to identify emission-line galaxies that are within the 
same dark-matter environment, and within a velocity range of 
— 1580 km s -1 to +3000 km s ' 1 from the HzRG. Typical galaxy 
overdensi ties around HzRGs ha ve velocity dispersions of <900 
km s _1 dVenemans et alJl2007 ). which suggests that the asym¬ 
metry of the velocity range in our case does not hinder our detec¬ 
tion of a potential overdensity of line emitters around the radio 
galaxy out to several times the velocity dispersion of the cluster. 
We present the results of both parts of this narrow-line imaging 
project in this section. 


5.1. Identification of candidate Ha emitters 


We identified candidate emission-line galaxies associated with 
the dark-matter environment of TXS 2353-003 candidate HAEs 
in two different ways. The first is a method introduced by 
[Bunker et al, ( 1995 ) that searches for a flux excess in the narrow 
compared to the broadband image, and the second is the direct 
inspection of the continuum-subtracted line image described in 


For the flux ex cess method of Bunker et al.1 ([1995), we used 
SExtractor v.2.5.0 (iBertin & Arnoutsl 19961) to construct catalogs 
from the narrow- and broadband images, setting the parame¬ 


ters DETECT_THRESH = 2.0 and ANALYSIS_THRESH = 2.0, while 


leaving all others at their default values. We then take the po¬ 
sitions of the galaxies in the catalog extracted from the broad¬ 


band image to perform aperture photometry within 2" aper¬ 
tures in both images. Using fixed apertures is appropriate for 
marginally or unresolved galaxies like ours, and has the advan¬ 
tage of providing errors that are independent of galaxy size (see 
also lKurk et al.l2004al) . Selecting our candidates from the broad¬ 
band image, rather than the narrowband, is different from many 
other narrowband searches, but makes our selection particularly 
robust, since it requires independent detections in two images. 
The disadvantage is that we might miss galaxies with particu¬ 
larly high emission-line equivalent widths. However, compari¬ 
son of the catalogs extracted from the two images shows that 
this was not the case in our analysis. 

We use the “significance of excess”, E, as defined by 
iBunker et al.1 ( 199 5). i.e., “the number of standard deviations 
between the counts measured in the broadband and the num¬ 
ber expected on the basis of the narrowband counts (assuming 
a flat spectrum)”, to select the besd Ha candidates. In addi¬ 
tion, and again following lBunker et al. ( 1995 ) we require that the 
sources fall above a given equivalent width threshold. Following 
iKurk et ahl (l2004al) . we choose EWo >50A and 25A. These 
thresholds correspond to broad- and narrowband colors of 0.38 
and 0.17, respectively. Three sources with E > 3.0 even have 
EWo > 100A, including the radio galaxy, which corresponds to 
a broad- to narrowband color of 0.69. 

The resulting color-magnitude diagram is shown in Fig. [9] 
We find seven sources with E > 2.0 and EWq > 50A. We also 
show another seven Ha candidates that have somewhat lower 
excess significances E =1.5-2.0. A fraction of these are probably 
part of the overdensity of TXS 2353-003. However, the fraction 
of misidentifications due to larger spectral slopes, for example, is 
likely to be larger among these sources, which is why we restrict 
our analysis to the most robust candidates above E =2.0. 


5.2. An overdensity of Ha candidates around TXS 2353— 003 

Given the relatively good seeing of our narrowband data of 0 A", 
most Ha candidates are spatially resolved (Fig. lA. 1 b and we es¬ 
timate Ha equivalent widths of 50—250 A in the rest frame, cor¬ 
responding to Ha fluxes between about 1 and a few x 10~ 16 erg 
s' 1 cm~ 2 . Their properties are listed in Table[2] We did not cor¬ 
rect the Ha flux and equivalent widths for [Nil] emission, which 
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Fig. 10. Positions of the candidate Ha emitters surrounding 
TXS 2353-003 on the sky. The HzRG is at the field center and 
labeled with the symbol “a”. Offsets are given relative to the po¬ 
sition of TXS 2353-003. North is up and east to the left. 


should fall into the same filter. Observations of UV/optically se¬ 
lected galaxies at high redshift typically show very low [Nil 1/Ha 
line ratios on the order of 10% (e.g.. (Forster Schreiber et afl 
12009 ?: Lehnert e t al .1 ,2009b lOuevrel et ali 120121) . Using the con¬ 
version of KennicuttT i 1998 ) with a 1 — 100 M 0 Salpeter initial 
mass function, we find star formation rates of 16-35 M 0 yr 1 per 
galaxy, and a moderate total star formation rate of 138 M 0 yr -1 
in these galaxies combined. These values were not corrected for 
extinction , so intrins i c star for mation rates could be higher. For 
example. iGarn et all (1201 ()l) and lStott et alJ (120 1 3h typically find 
1 mag extinction in H-alpha meaning that intrinsic star formation 
rates are probably a factor of ~2.5 times higher. The expected 
fraction of low-redshift interlopers is much lower in our case 
than for narrowband searches using Lya. The brightest emission 
lines longward of Ha are the Paschen and Brackett hydrogen 
lines in the near-infrared, with typically much lower equivalent 
widths than Ha. 

To quantify whether the observed number of Ha candidates 
around TXS 2353-003 may herald an overdensity of galaxies, 
we have to compare our data with samples of Ha emitters in 
the field at the same epoch. A blind redshift search for HAEs 
at z = 1.48 has recently been performed bv ISobral et al.1 ( 2013 ?) 
through the NB 1.617 filter at UKIRT, which has a very sim¬ 
ilar width of 210 A (compared to 250 A for our ISAAC fil¬ 
ter), that we use as reference f or the expected source density 
in the field. The HiZEL survey dSobral etal .1120091120 1 ()L 120121 
1201 3l) gives an estimate of the Ha emitters that we can expect 
in a blank field: ISobral et al.l ( 2012 1 concentrated on Ha emit¬ 
ters at z ~ 1.47, very close to the z=1.49 of TXS 2353—003. 
They found 295 NB emitters over the 0.67 deg 2 field of their 
main analysis, and 411 over the entire 0.79 deg 2 field of their 
NB observations. Among them, 190 objects also show an ex¬ 
cess in deep NB imaging centered on the [Oil]A3727 emis¬ 
sion line and have colors compatible with a photometric red¬ 


shift z ~ 1.47. This leads to densities of candidate HAEs in the 
field of Lhae — 0.14 HAE arc min 2 (with all NB emitters) ; 
T-hae — 0.12 HAE arc min 2 (with NB emitters in their main 
field) and 'Lhae — 0.08 HAE arc min 2 (with NB emitters having 
an adequate photometric redshift). These values include galaxies 
down to a 3a flux level of 7 x 10 -17 erg s -1 cm -2 , compared to 
1 x 10 -17 erg s -1 cm* 2 in our case. 

Around TXS 2353—003, with an effective field of view of 
5.4 arcmin 2 (neglecting the trimmed edges where our data did 
not reach the full depth) we find six Ha emitting candidates with 
fluxes > 7 x 10 -17 erg s -1 cm* 2 and significances E > 2.0 in 
addition to the radio galaxy. This corresponds to a source density 
of 1.1 arcmin -2 , and an overdensity by a factor ~8 relative to the 
field. 

iKurk et aki (2004bh detected 28 HAEs in two ISAAC point¬ 
ings around MRC 1138-262 at z = 2.16 down to an Ha flux limit 
of 2.5 x 10 17 erg s * 1 cm -2 , which corresponds to the same lu¬ 
minosity limit as the current analysis. They found a source den¬ 
sity of 2.2 arcmin -2 . Although our source density at z=1.5 is 
lower than that in the z=2.16 overdensity around MRC 1138- 
262 in absolute terms, the significance of the overdensity relative 
to the field at the same epoch is high, because the star formation 
rate density at z = 1.5 in the field is already stron gly declining 
for the highest Ha luminosities (e.g.. ISobral et al.112013i> . 

Finding high surface-brightness gas out to the periphery of 
our SINFONI data cube raises the question whether there is sim¬ 
ilarly bright line emission outside the small, 8" x8” field of view 
of SINFONI (67 kpcx67 kpc at z=1.49). We do not find any 
such emission down to a 3a detection limit of 6 x 10 20 erg s 1 
cm -2 A -1 in the continuum subtracted narrowband image. With 
the 250A width of the NB1.64 filter, and cosmological surface- 
brightness dimming by a factor (1 +z) 4 = 39, we are sensitive to 
the equivalent of > 60 x 10 -17 erg s -1 cm -2 arcsec -2 in nearby 
clusters. This is more than an order of magnitude fainter than the 
surface brightness observed, e.g., in parts of the filaments in the 
Perseus cluster (lHatch et alJ(2007h . 


6. Comparison of TXS 2353-003 and 

NVSS J210620-214003 with cluster central 
galaxies 


High-redshift radio galaxies do not often reside in solitude, 
but are surrounded by overdensities of line or continuum emit¬ 
ters, suggesting that a significant fraction of the HzRG popu¬ 
lation may be the proge nitors of the central galaxies of mas¬ 
sive galaxy clusters (e.g . . [Chambers et al. l T99GL LeFevre^tal 


1996; Kurk et al. 2004a Venemans et a 1 .1 l2007b G a I a me tz et al 


2012). Wvlezalek e t al] (12013 ) obtained Spitzer warm-mission 
IRAC photometry at 3.6/im and 4.5um of 200 HzRGs, includ¬ 
ing our two sources. They find thafi_ compared to the Spitzer 
UKDSS Ultra Deep Survey dKim et al.ll201~fl) . 55% of HzRGs 
are surrounded by > 2a overdensities of galaxies with colors 
between the 3.6um and 4.5/im channels that are consistent with 
redshifts z>1.3. Like any method identifying galaxy overdensi¬ 
ties, however, the IRAC color selection only samples parts of the 
population of putative satellite galaxies around HzRGs, and can 
therefore only provide lower limits to the actual density of such 
galaxies around a particular HzRG. It can therefore only demon¬ 
strate the presence, but not the absence of such an overdensity of 
gal axies around a gi v en r adio galaxy. 

IWvlezalek et d] (hoi.ll) show that TXS 2353-003 is sur¬ 
rounded by one of the highest overdensities of IRAC-selected 
sources, with an excess of > 5a, and consistent with the re- 
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suits of our narrowband search presented in ij5] NVSS J210626- 
314003, however, does not stand out as a source with a par¬ 
ticularly dense environment of IRAC-selected sources in the 
tWvlezalek et al.l ( 2013 ) sample, a 1.2 ff excess. Unfortunately, 
the observing program committee of ESO only granted us nar¬ 
rowband imaging^of TXS 2353-003, not of NVSS J210626- 
314003, so that we were unable to do a narrowband search 
around this galaxy. We note, however, that selecting a galaxy 
based on its bright radio emission already leads to sign ificant 
biases toward galaxies in dense environments ( e.g., iBes 2 [2000:; 
Ramos Almei da et al.l 120131 lHatch et al.l 120141) . and we argue 
in the following that the stellar component of NVSS J210626- 
314003 itself shows some of the characteristic signatures of clus¬ 
ter central galaxies. 


6.1. K—z relationship 


iKauffmann & Chariotl (il998l) argued that the observed K-band 
magnitude can be used to approximate the stellar mass of all but 
the most rapidly growing galaxies out to redshifts z>2, with a 
scatter of about a factor of 2. This is consistent with the empiri¬ 
cal tight K-z relationship between observed K-band magnitudes 


of powerful radio galaxies and their re dshifts (Lilly & Longair 
1984 iDe Breuck etafl 120021: IWillott et al.l 120031: iBrvant etal 


20091) and the small mass range of these galaxies in the Spitzer 
photometric survey of HzRGs by ISevmour et al.l ([2007) and 
IDe Breuck et al.l (l2010l) . 

In Fig. E] we show where TXS 2353-003 and 
NVSS J210626—314003 fall re lative to the K - z rela tion- 
( 20071) and IDe Breuck et al.l (l2010h . and 
20121) , who found a similar relationship 


Seymour et al 


Lidman et al. 


ship of 
that of 

for brightest cluster galaxies (BCGs). This figure shows the 
good agreement of our targets with both samples, as well as the 
self-consistency of both samples with each other (as expected if 
HzRGs are the progenitors of BCGs). 

Encouraged by th e re sults of [K auffman n & Chariotl (JT998); 
ISevmour et al.l d2007l) , and lDe Breuck et al. ( 2010I) . which, when 
considered together, suggest that the observed K-band magni¬ 
tude of HzRGs should scale with the stellar mass, we used 
the K-band magnitudes of z ~ 2 HzRGs in the ISevmour et al.l 
(120071) sample and the stellar mass estimates given in the 
same paper, to obtain a simple empirical relationship between 
stellar mass, M sle u ar , and observed K-band magnitude, m.K, 
\og{M ste ii ar /M q) = 6.18 — 0.189 nix . We find ala scatter of 
about 0.3 dex in mass, consistent with that previously found 
bv lKauffmann & Chariotl ( 1998 ). We note that we are using ob¬ 
served magnitudes from galaxies in a similar redshift range, and 
therefore do not need a k-correction. 


This allows us to derive a rough mass estimate for our two 
sources. For TXS 2353-003 and NVSS J210626-314003, K- 
band magnitudes of 18.1±0.3 mag and 18.7±0.2 mag suggest 
stellar masses in the range log (M 2353 /M e ) = 11.1 — 11.4 and 
log (M 21O6 /M 0 ) = 11.2—11.5, respectively. While such an ap¬ 
proach should not replace more detailed analyses of the con¬ 
tinuum emission where additional photometric constraints are 
available, it does highlight that our two galaxies fall into the typ¬ 
ical range of stellar masses of powerful radio galaxies at these 
redshifts. 


6.2. Surface-brightness profiles 

In Section 14.21 we showed that both galaxies are spatially re¬ 
solved in our ground-based imaging, and that both have surface- 



Fig.ll . Observed K —band magnitude of brightest-cluster galax¬ 
ies (shown as '+’) and radio galaxies (shown as ‘ x ’) as 
a function of redshift (plot reproduced from iLidman et ali 
I2012L and slightly extrapol ated to z=2.5). S ome of the data 
points were taken from IStott et all d2008l) and IStott et all 
(120101) . The large red stars illustrate that TXS 2353—003 and 
NVSS J210626—314003 fall well within the expected range of 
K —band magnitude, and also have typical K-band magn itudes 
for HzR Gs at their redshifts (black crosses Seymour et al .1120071 
see also IDe Breuck et al.ll2002t IWillott et alJl2003 ; Brvant et al.l 
2009). At z> 2, Ha and [NII]AA6548,6583 fall into the K-band 


filter and can affect the continuum magnitudes by up to about 
0.3 dex in extreme cases (iNesvadba et alJi2()()6l) . 


brightness profiles that are not consistent with the pure Sersic 
profiles of early-type galaxies (n=4). The central arcsec of 
the profile of NVSS J210626—314003 can be fitted with a de 
Vaucouleurs profile with R=15 kpc, but fainter extended emis¬ 
sion persists at larger radii. TXS 2353-003 can only be fitted 
with a de Vaucouleurs profile when assuming an unphysically 
large radius of 26 1 kpc . 

iTargett et alJ (1201 It) analyzed the K-band morphologies of 
13 radio galaxies at z=1.5—2.5 observed with UKIRT during 
very good seeing, and found that most sources are well fit with 
Sersic indices n= 3—5 and effectiv e radii between 5 and 10 kpc. 
This is in contrast to the findings of lPentericci et alJ ( 2001 ) with 


HST/NICMOS through the F160W filter that only 5 out of 19 
sources have well-established de Vaucouleurs profiles (although 
the iPentericci et al.l 12001). sam ple does extend to higher red¬ 
shifts than that of Targett et al.l 1201 ll and, with H-band obser- 
vation s, samp les much shorter wavelengths). Only one source of 
the iTargett et al.l (1201 ll) sample, 0128—264, has an effective ra¬ 
dius R=15 kpc, as does NVSS J210626-314003. This galaxy is 
also included in our parent sample with SINFONI data. It has 
a small emission-line region, although elongated along the radio 
jet axis, and a very extended (294 kpc) radio size, resembling the 
two sources discussed here (Nesvadba et al. 2015, in prep.). 

Massive early-type galaxies at high redshift appear to be 
generally characterized by their greater compactness com- 
pared to similar, more nearby galaxies dDaddi et al.l {2005; 
Ivan Dokkum et all 120081 who found R e < 1 kpc and 0.9 kpc, 
respectively ). Typically, these galaxies are only resolved with 
deep, high-resoluti on HST imaging, showing no deviations from 
Sersic law s (e.gJSzomor u et all COl OHCassata et alJl2010l) . The 
HzRGs of Ta rgett et al.l (1201 ll) . which have relatively large radii 
(average 8 kpc) compared to submillimeter galaxies (average 
3 kpc) at similar redshifts, show an average decrease in size by 
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a factor 1.5 compared to equally massive galaxies at low red- 
shift. A size decrease by a factor of ~ 2 has also been found by 
iPentericci et al.l (1200 ll ) for their HzRGs that are well modeled 
with a de Vaucouleurs profile. 

Extended stellar envelopes are a primary characteristic of 
central cluster galaxies in the low-redshift Universe (cD galax¬ 
ies; e.g.. lSchombertlll987l) . and likely formed from the debris 
of repeated accretion of several satellite galaxies. This inter¬ 
pretation also agrees with the finding of Targett et alJ (1201 ll) 
that HzRGs fall onto the same Kormendy relationship between 
central surface brightness and effective radius as massive low- 
redshift galaxies, which, as stated by the authors, would imply 
that the size increase from z~2 to today must be accommodated 
by an increase in stellar mass at large galactocentric radii, as 
could be produced by repeated accretion of satellites. iBest et all 
( 1998 ) identified similar structures in a stack of 12 radio galax¬ 
ies from the 3CRR sample at z= 0.6 — 1.4, but to our knowledge, 
the only example of a HzRG with extended stellar halo at greater 
redshift currently known in the literature is MR C 1138—262 
at z=2. 2 (the Spider web Galaxy. iMilev et akll2006h . lHatch et alJ 
(2008|; lHatch et aD (120091) discussed the envelope around this 
source in detail, finding that the rest-frame UV colors of the faint 
continuum emission are consistent with a stellar envelope (rather 
than scattering from a dusty gaseous halo) with a formation his¬ 
tory that is coeval with that of the high-surface brightness central 
regions. This is in contrast to semi-analytical simulations which 
postulate that such envelopes assembled gradually through dry 
merg ers over cosmological timescales (e.g. JDe Lucia & Blaizofl 
120071) . 


7. The nature of the extended ionized gas 

A good alignment between jet and g as is ty pical of distant radio 
galaxies (the alignment effect; e.g., Cimatti et al.lll993|). and this 
alignm ent was also among the prime ar guments of [Collet et akl 
(1 201 5l) and of INesvadba et alJ (120061) : Nesvadba etal] (120081) 
that jet and gas are physically related in their SINFONI samples 
of HzRG. They also found very large, super-gravitational veloc¬ 
ity gradients reaching to more than 1000 km s 1 associated with 
the most powerful radio sources, and broad line widths in the 
extended gas with FWHM > 800 km s' 1 . They considered this 
as evidence of fast, AGN-driven outflows in HzRGs. 

None of this applies here. Velocity gradients are well below 
1000 km s , line widths outside of the central regions are only 
on the order of 200-300 km s , and the mismatch in position 
angle between jet and gas is striking. This mismatch makes it dif¬ 
ficult to postulate that jet cocoons are driving the gas out as they 
expand through the ambient gas, because it appears unavoidable 
in this model that cocoon (and hence, jet) and gas are cospatial, 
both in three dimensions and in projection. 

The presence of this off-axis component of line emission is 
worth discussing, as is the absence of bright, extended emission¬ 
line regions along the jet axis. TXS 2353—003 may have a faint, 
not well resolved emission-line component along the jet axis and 
within the galaxy itself. However, this does not stand out com¬ 
pared to the gas in other regions of the galaxy, and is therefore 
difficult to associate uniquely with the radio jet. We find no sim¬ 
ilar emission at larger radii, including the cluster scales that we 
cover with ISAAC around TXS 2353-003 (©. 

The absence of bright gas on scales much larger than the 
host galaxy cannot be a mere observational effect, since the sen¬ 
sitivities of these data are comparable to the overall sample, and 
the surface brightness of the line emission in our targets is not 
significantly lower than in other HzRGs with SINFONI data. 


However, for the largest sources with ’regular’ line emission, in 
particular MRC 2104-242 at z=2.5, the extended line emission 
can become very faint and is concentrated along long, thin, very 
filamentary structures, suggesting that the warm gas is either be¬ 
ing heated to temperatures T3> 10 4 K, or becoming strongly dis¬ 
persed as the cocoon continues to expand (Nesvadba et al. 2015, 
in prep.; see also iPentericci et akl1200 ll who observed the same 
gas as faint, extended filaments in their HST/NICMOS broad¬ 
band imaging). 

We stress that finding extended gas around HzRGs that 
is unrelated to the radio jets is by itself not uncommon. 
Many HzRGs are surrounded b y diffuse halos of warm ionized 
gas (IVillar-Martln et akl 120031) . however, dense, high-surface- 
brightness Lya emission consistent with the Ha surface bright¬ 
nesses that we o bserve here is only found w ithin the jet cocoon 
in the galaxies of lVillar-Martm et ak ( 2003 ) Ha associated with 
the diffuse ionized gas would be about 2 orders of magnitude 
fainter than what we observe. 

We will now discuss several hypotheses, from AGN orien¬ 
tation to galaxy interactions to extended gas disks or filaments 
within the radio galaxy itself or its immediate surroundings, to 
explain the nature of this gas. 

7.1. Partial AGN illumination of ambient halo clouds 

The line ratios of both sources suggest that they are photoion- 
ized by the central AGN (fj3). This could imply that they do 
not represent a distinct structure, but are merely part of a 
population of ambient clouds that are distributed over large 
solid angles, and of which a subset is being lit up as it inter¬ 
cepts the quasar ionization cone. Filaments of neutral gas have 


also been observed around some HzRGs (Ivan Oiik et al. 

1997 

Jarvis et al. 2003|Wi 

Iman et al. 2004; De Breuck et al. 

2005 

Nesvadba et al.E009; 

Emonts et al. 20141). 


However, not all HzRGs _show signat ures of neutral gas 
clouds in their halos. Ivan Oiik et ak ( 1997 ) found Lya absorp¬ 
tion only in galaxies with radio sizes <50 kpc, in stark con¬ 
trast to the large radio s izes, 200 kp c, of TXS 2353-003 and 
NVSS 210626-314003. iBinette etakf ( 2000 ) showed from an 
analysis of the ionization properties that these absobers must 
also lie at larger radii than the radio source. In the more com¬ 
pact radio galaxies, Lya absorbers are found not only against the 
nucleus, but also_against more extended Lya emission around 
the galaxy (iKurk et al.1 l2004al) . which suggests that this result 
holds for large solid angles around the radio galaxies, and not 
just along the radio jet axis. 

This scenario would also require a misalignment between 
radio jet axis and quasar illumination cone, as is possible be¬ 
cause of jet precession, or a misalignment between jet axis and 
the normal of the torus, for example. iDrouart et ali (l2014i) re¬ 
cently argued, based on Spitzer mid-infrared imaging and the 
core-to-lobe fraction of the radio sources, that the unified model 
must hold for the global population of HzRGs, but the presence 
of nuclear broad Ha emission in a few sources does also sug¬ 
gest that this might no t be the case for each individual HzRG 
(INesvadba et al.ll201 lal) . 

Alternatively, we may be seeing gas in an old cocoon from 
a previous radio-loud activity phase of the super-massive black 
hole. Jet precession has been invoked to explain the X-shaped 
radio sources at low redshift, where about 5-1 0% o f FRII radio 
sources have two pairs of radio lobes dLeahv & Muxlowll2002l) . 
The second bar of the “X”, typically secondary pairs of fairly 
diffuse, low surface-brightness wings in centimeter radio con¬ 
tinuum imaging, may be the relic of a previous burst of nuclear 
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activity along a different jet axis. Detecting such low surface- 
brightness emission directly would be very challenging at high 
redshift. 

This scenario, however, also has its difficulties. 
iKaiser & Cotter! ( 2002 1 argued that clouds in relic cocoons 
will be destroyed by shocks within a few 10 6 yrs, which is of the 
same order a s the ages of bright radio so urces in high-redshift 
galaxies (e.g jBlundell & Rawlings! 19991) . In the cocoon model, 
dense clouds are confined by the high pressure of the cocoon 
material, and should disperse within a sound-crossing time 
once the jet has stopped maintaining the cocoon pressure high 
(e.g., iFabian et al.l 119871) . From the size of the two lobes in 
our sources, 310 kpc and 190 kpc, and assuming a jet advance 
speed on the order of 0.1c, we can roughly constrain that the 
feeding of a putative older jet component must have ended at 
least about 1 x 10 7 yrs ago in our two sources, so that most of 
the emission-line clouds should have already been destroyed 
or evaporated. The relatively old age of our two sources does 
not make them good candidates for seeing relic cocoons, in 
particular since many younger radio sources, which should in 
principle have brighter relic cocoons, do not show evidence of 
two sets of radio lobes or extended emission-line regions. 

We conclude that neither partial AGN illumination of a gen¬ 
eral population of ambient clouds nor repeated cycles of radio- 
loud AGN activity with different jet orientations appears to be a 
good explanation for the nature of our sources. 

7.2. External gas supply from a satellite galaxy? 

Gas transfer from a satellite galaxy undergoing accretion onto 
the HzRG, a s prev iously observed in the z=3.8 HzRG 4C60.07 
(ilvison et al .1120081 ). is another interesting hypothesis for the ori¬ 
gin of this misaligned gas in our two sources. In this case, the 
extended line emission could trace a gas tail produced by ram- 
pressure stripping or tidal forces. However, only TXS 2353—003 
has a companion along the direction of the extended gas, but at 
a larger distance from the radio galaxy. NVSS J210626-314003 
has no obvious companion within several tens of kpc that can be 
associated with the extended emission-line region. 

It is uncertain whether the galaxy near TXS 2353-003 
(Fig. 1) is physically associated with the radio galaxy, although 
its infrared colors as measured with IRAC are consistent with 
a redshift z<;1.3 using the criteria of IPanovichl (120081) ((0). 
However, we detect no Ha line emission down to 3.3 x 10~ 18 
erg s 1 cm~ 2 . Tidal forces during galaxy interactions produce 
extended tails, and also funnel gas t oward the nuclei of the galax¬ 
ies (e.g.. [Barnes & Hernuuisti 1996i) , fueling intense nuclear star- 
bursts. Likewise ram pressure should not only strip parts of the 
ISM, but also compress the gas along the head of the infalling 
galaxy, enhancing the gas surface brightness within the galaxy, 
and potentially star formation dKanferer et a1.1(20081) . 

It appears therefore difficult for a galaxy to produce a tail 
of ~ 10 9 M 0 of ionized gas without sustaining significant star 
formation. Our upper limit of Ha flux corresponds to star for¬ 
mation rates on the order of 10 M 0 yr 1 (using the iKennicuttl 
119981 . calibration and neglecting extinction). Since the source is 
bright in the continuum at similar wavelengths, the Ha equiva¬ 
lent width must also be low. Although gas, dust, and young stars 
are not necessarily co-spatial dCalzettill 19971) . this would imply 
here that the galaxy, in spite of losing considerable fractions of 
its ISM in the putative tail, is also maintaining a highly efficient 
obscuring dust screen around its star-forming regions. This ap¬ 
pears contradictory, in particular for high-z galaxies, which typ¬ 
ically have extended star formation in multiple knots across the 


galaxy (e.g.. [Forster Schreiber et aD 12009 ). Other HzRGs, e.g., 
MRC 1138-262 at z = 2.16, are already known to be surrounded 
by companions within a_few tens of kpc that are already on 
the red sequen ce (Koda ma et al.l 120071) an d lack line emission 
dNesvadba et al.l l2006h iKuiner et al.l 1201 ll) . It thus appears un¬ 
likely that external gas supply from a satellite galaxy is a good 
explanation of the extended emission-line regions in our galax¬ 


ies. 


7.3. Extended gas disks within the radio galaxy 


iBesd ( 2000 ) pointed out that radio galaxies with very extended 
radio sources are particularly good candidates for being clus¬ 
ter central galaxies because the pressure from the Mpc-scale 
intracluster medium boosts the luminosity of the radio source 
even for comparabl y extended jets (see also lAthreva et al.l l 1998: 
Kla mer et al.l 12006). Many authors have previously suggested 
that the progenitors of brigh test cluster galaxies are likely to 
be found among HzRGs (e.g., 

12006; I lalch et al.ll2009l) . and 


Pentericci et al 


Wvlezalek et al 


2001 1 Milev et al. 


(2013) find from 


IRAC imaging that TXS 2353-003 is even associated with the 
most strongly pronounced overdensity within their sample. 

The gas that we see here may actually have much in com¬ 
mon with the extended gas disks or filaments that are found in 
and near cluster central galaxies in the more nearby Universe. 
Mismatches between the position angles of jet and gas are fairly 
common in about 30% (iMcDonald et al.l 120108) of cool-core 
galaxy clusters with extended Ha filaments. A particularly clear 
example, where the warm ionized gas appears to avoid the cavi¬ 
ties that the radio jet has inflated in the intracluster medium is 
Abell 1795, in particular_at the highest emission-line surface 


brightnesses (Ivan Breugel et alJ 1984 ). In total, su ch structures 


have up to a few 10 10 M 0 (e.g .. [Salome et alfeOObl) of warm and 
cold gas with T< 10 4 K, typically dominated by cold molecular 
gas. We have already argued in lj6] that our sources follow the 
overall trends found in BCGs between Ha luminosity and sur¬ 
face luminosity, size, and ratios of low-ionization emission lines. 


tical line emission (e.g.. 

Hatch et al. 2006; Wilman et al. 

2006; 

Hatch et alJl2007: Wilman et al. 

2009; Farage et alJl2010L 

20121) 


shows that the gas has typically FWHM ~ 100 — 300 km s' 1 
outside the near-nuclear regions, with a characteristic broaden¬ 
ing toward the center, with typical FWHM > 500 km s 1 and up 
to the values we find near the center of NVSS J210626—314003 
and TXS 2353-003 (0 and Figs. 0 and 0). However, reg¬ 
ular velocity gradients of a few 100 km s -1 , similar to 
what we observe in our sources, are found only in some 
BCGs, e.g., Abell 262 dHatch et alJ 120071) or 2A 0335+096 
dFarage et al.ll2() 12 ). Alternatively, disks with well-ordered ve¬ 
locity fields, perhaps more like in NVSS J210626-3140003 
than in TXS_2353—003, have also been found, for example in 
Hydra A (iHamer et alJ 12014 ). Such disks are also typically not 
aligned with the radio jet axis, but can be illuminated by the 
central AGN, in particular if the opening angle of the AGN illu¬ 
mination cone is large. 

The observed velocity gradients in our two sources are in 
fact consistent with rotational motion in a gravitational potential 
as suggested by the stellar mass estimates of 1.5 — 3 x 1O U M 0 
(ED- A dynamical mass estimate, v, can be derived by setting 
v = \JM sin i-G /R. where M is the stellar mass, R the radius of 
the putative disk, i the inclination angle, and G the gravitational 
constant. With R=10.5 kpc, measured for the very regular veloc¬ 
ity field of NVSS J210626-314003, we find v= 250 - 350 sin i ~ 1 
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km s 1 , compared to 220 km s 1 observed. Even if the gradi¬ 
ents were dominated by radial outflow motion, the small velocity 
range implies that most of this gas is unlikely to escape. 

This is consistent with the finding that the gas extends over 
similar radii to the faint continuum wings in the ISAAC im¬ 
age of TXS 2353-003 on both sides, and in NVSS J210626- 
314003 on the southwestern side. If the wings in the contin¬ 
uum surface brightness profiles of our two sources represent 
extended stellar halos (lHatch et al.l 120091) . perhaps originating 
from a phase of rapid accretion of man y sat ellite galaxies in the 
very early evolution of the HzRGs dBurkert et al.ll2008l) , then we 
may be seeing leftover gas that is settling down after this phase. 
[McDonald et alJ ( 2010 1 suggested that the condensation of such 
gas into fairly dense clouds could be accelerated by weak shocks 
caused by the passage of a companion galaxy, as may be the case 
for the putative companion in TXS2353-003. 

Weak shocks could also be produced by the vestiges of 
the relic cocoon that initially formed when the jet was passing 
through the inner regions of the halo. This scenario has some 
resemb lance to the cyclical AGN f eedback scenarios invoked, 
e.g,, b y lPizzolato & Sokerl (l2005h or lAntonuccio-Delogu & SiTkl 
(1201 fll) to explain exten ded gas disks and filaments in nearby 
cluster central galaxies. iNesv adba et al] (1201 Ol 1201 lbl) argued 
that isolated radio galaxies with extended gas disks may also 
require repeated, e pisod ic AGN activity, broadly akin to the 
[ Pizzolato & Sokerl (120051) model, in order to understand the ob¬ 
served properties of the gas. Winds associated with the large ra¬ 
diojets we observe may have dispersed the ISM of these galax¬ 
ies just a few 10 7 yrs ago, and parts of this gas may now be 
raining back onto the galaxy, possibly feeding a new feedback 
cycle. The age of the Universe at z= 1.5-2 is already > 3 Gyr, 
so that repeated activity cycles are possible for duty cycles 
of a few times 10 8 yrs, as seems appropriate at low redshift 
(iPizzolato & Sokeri[2005l) . The gas cooling times could be low¬ 
ered by compression of the outflow and diffuse ambient gas that 
is being swept up as the two bubbles of the AGN cocoon con¬ 
tinue to inflate in the form of a mo mentum-driven wind af ter the 
feeding from the AGN has ceased dKaiser & Cotterli2002l) . This 
cyclical feedback may contribute to forming the hot intraclus¬ 
ter X-ray gas in massive galaxy clusters today, which presum¬ 
ably too^_place_around_redshift z>2 dNath & Rovchowdhurvl 
12 0021; McCarthy et afl2008i) . Energy injection through repeated 
episodic AGN activity has also been invoked as an explana¬ 
tion of how the X-ray halos surrounding individual massive 
early-type galaxies can be maintained over a Hubble time (e.g., 
iMathews & Brighentill2003liBest et aDl2006h . 

8. Summary 

We have presented an analysis of rest-frame optical imaging 
spectroscopy and deep broadband near-infrared imaging of two 
radio galaxies at z~2, which have extended emission-line re¬ 
gions that are strongly misaligned relative to the axis of the ra¬ 
diojets. This is in stark contrast to the majority of powerful radio 
galaxies at similar redshifts, where the gas is aligned within 20- 
30° of the jet axis, as expected from a cocoon of turbulent gas 
that is being inflated by the expanding radio jets. For one source, 
TXS 2353-003 at z=1.5, we also present Ha narrowband imag¬ 
ing through the [Fell] 1.64 filter of ISAAC, finding an overden¬ 
sity of Ha candidates which resembles those of radio galaxies at 
higher redshift. 

The gas in both galaxies is less perturbed than in radio galax¬ 
ies with more typical gas properties. In particular, line widths in 
the extended gas are lower and, with FWHM ~200-300 km s 1 


comparable to those in the extended gas disks and filaments sur¬ 
rounding brightest cluster galaxies at lower redshifts. Overall, 
we find remarkable similarities to the BCGs in low-redshift cool- 
core clusters, ranging from the extended, faint continuum halos 
in both galaxies that are atypical for massive high-redshift galax¬ 
ies except for the particularly well-studied Spiderweb Galaxy 
MRC 1138-262, where such structures are interpreted as ex¬ 
tended stellar halos broadly akin to cD galaxies at lower red¬ 
shifts. Both galaxies fall onto the relationship between K-band 
magnitude and redshift for brightest cluster galaxies and HzRGs, 
suggesting they are among the most massive galaxies at their 
epoch. We note that, given that HzRGs are now generally con¬ 
sidered to be the progenitors of BCGs, the distinction between 
both classes is not mutually exclusive, and our two sources may 
simply be further evolved than many other HzRGs, but still fol¬ 
lowing the same overall evolutionary path. 

We discuss several scenarios for the nature of this gas, find¬ 
ing that simple illumination effects in otherwise typical HzRG 
halos is not a good explanation, and neither are galaxy inter¬ 
actions. Generally speaking, our results support the hypothesis 
that the extended line emission arises from extended gas disks 
or filaments within or near the radio galaxy, by analogy with 
broadly similar structures in nearby cluster central galaxies. This 
could be evidence for cyclical AGN feedback, which has been 
discussed on several occasions for nearby clusters. 
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Table 1 . Properties of the two galaxies discussed in this study. 


TXS 2353-003 NVSS J210626-314003 

Redshift 

1.487 

2.104 

Dz, [Gpc] 

10.8 

16.6 

radio power [10 28 W Hz~’] 

3.7 

4.2 

radio size [arcsec] 

38.8 

24.2 

radio size [kpc] 

328 

201 

K-band magnitude [mag] 

18.1T0.3 

18.7T0.2 

PA offset between jet and gas 

~ 90° 

~ 60° 

F([OIII]A5007) fl 


66± 13 

F(Hj3) fl 


5.2± 1 

F([NII]A6548f 

2.8±0.6 

1.2±0.2 

F([NII]A6583) fl 

8.3 ± 1.7 

3.7 ±0.7 

F([SII]A6716) fl b 

1.4±0.5 

2.3±0.5 

F([SII]A6731f 

4.1 ± 1.0 

2.2±0.5 

F(Ha)° 

12.0± 2.4 

14.7 ± 3.0 

±?(Ha) [10 43 ergs^'] 

1.7 

4.8 

ionized gas mass [ 10 8 M 0 ] 

2.4 

9.3 

SINFONI rms c 

4.3 

6.2 


All fluxes are given in units of 10 16 erg s 1 cm 2 . W [S1I]A6716 in TXS 2353—003 is affected by a telluric line. Surface brightness 
detection limit in units of 10~ 17 erg s~’ cm~ 2 arcsec~ 2 . 


Table 2. Parameters of the putative Ha emitters around TXS 2353-003 found in our observations. 


ID 

RA 

(J2000) 

Dec 

(J2000) 

distance 

[arcmin] 

NB flux a 
[O.IpJyl 

BB flux a 

[AtJy] 

BB-NB color 

EW 0 b 

[A] 

Ha flux' 

[10 17 erg s 1 cm- 2 ] 

SFR 1 ' 

[M@ yr '] 

a 

23:55:35.58 

-00:02:46.2 

0.0 

25.9 

14.3 

0.69 

152.4 

63.6 

N/A 

b 

23:55:37.95 

-00:01:36.8 

1.29 

16.2 

11.3 

0.95 

100.2 

32.1 

35 

c 

23:55:32.56 

-00:01:59.2 

1.08 

10.3 

6.2 

0.86 

132.3 

23.8 

26 

d 

23:55:37.64 

-00:03:23.9 

0.81 

8.3 

5.4 

0.76 

112.7 

18.7 

20 

e 

23:55:37.74 

-00:01:35.3 

1.29 

6.2 

3.5 

0.91 

145.1 

20.1 

22 

f 

23:55:31.92 

-00:01:33.8 

1.51 

6.7 

3.5 

1.01 

166.8 

17.1 

19 

g 

23:55:30.72 

-00:01:45.7 

1.57 

6.8 

2.7 

1.3 

248.2 

15.0 

16 


Measured in 2" diameter apertures. Rest-frame equivalent width. ^ Assuming that the color excess in the NB image is entirely from line 
emission. ^ Following lKennicuttl d 1998l f. We do not give a star formation rate (SFR) for the radio galaxy, because the line ratios are not consistent 
with photoionization from young stars. 
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Fig.A.2. Ha emitters around TXS 2353-003 selected by the Bunker plot and confirmed when subtracting the continuum from 
the NB image. Left column: //—band images of the candidates. Middle column: NB1.64 images of the candidates. Right columns: 
Continuum-subtracted NB images of the candidates, the last column showing them smoothed with a Gaussian kernel of a = 1 pixel. 
The color palette spans —3(7 to +5(7 in each case in order to highlight the faint features of these HAE candidates. The contours 
give the morphology of the sources at larger surface brightnesses: white contours are at levels 3 ; 5 ; 10 ; 15 ; 20 <7 and the black 
contours are at levels —10 ; —5 ; —3(7. The HzRG is labeled a. 
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